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We have recently seen new upper bounds for B® — > a key decay 

to search for physics beyond the Standard Model. Furthermore a non- 
vanishing decay width difference Ar s of the B s system has been measured. 
We show that Ar s affects the extraction of the — > fi + fi~ branching ratio 
and the resulting constraints on the New Physics parameter space, and give 
formulae for including this effect. Moreover, we point out that AT S provides 
a new observable, the effective B® — > fi + fi~ lifetime t m + m -, which offers a 
rN theoretically clean probe for New Physics searches that is complementary 

to the branching ratio. Should the B® — > branching ratio agree with 

the Standard Model, the measurement of t^+ m -, which appears feasible at 
upgrades of the LHC experiments, may still reveal large New Physics effects. 
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We have recently seen new upper bounds for B° — > fi + /i~ , a key decay to search for physics beyond 
the Standard Model. Furthermore a non- vanishing decay width difference AF S of the B s system has 
been measured. We show that Ar s affects the extraction of the — > fi + /i~ branching ratio and 
the resulting constraints on the New Physics parameter space, and give formulae for including this 
effect. Moreover, we point out that AF S provides a new observable, the effective B° — > lifetime 
r^+^- , which offers a theoretically clean probe for New Physics searches that is complementary to 
the branching ratio. Should the B® — > fi + /i~ branching ratio agree with the Standard Model, the 
measurement of t m + m - , which appears feasible at upgrades of the LHC experiments, may still reveal 
large New Physics effects. 
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I. INTRODUCTION 

Thanks to the Large Hadron Collider (LHC) at CERN 
we have entered a new era of particle physics. One of the 
most promising processes for probing the quark-flavor 
sector of the Standard Model (SM) is the rare decay 
£?S? — > In the SM, it originates only from box 

and penguin topologies, and the CP-averaged branching 
ratio is predicted to be [I] 

BR(B S -> /j+aOsm = (3-2 ± 0.2) x 1(T 9 , (1) 

where the error is fully dominated by non-perturbative 
QCD effects determined through lattice studies. The 
most stringent experimental upper bound on this branch- 
ing ratio is given by BR(_B S — > fj, + n~) < 4.5 x 10~ 9 at 
the 95% confidence level (C.L.) 0. 

In the presence of "New Physics" (NP), there may 
be additional contributions through new particles in the 
loops or new contributions at the tree level, which are 
forbidden in the SM (see Ref. pQ and references therein). 

A key feature of the i3 s -meson system is £?°-£?° mix- 
ing. This quantum mechanical effect gives rise to time- 
dependent oscillations between the B° and B® states. In 
contrast to the Bd system, we expect a sizable difference 

(s) (s) 

Ar s = — Th between the decay widths of the light 
and heavy B s mass eigenstates [3]. 

Performing a time-dependent analysis of B Q S — > J/ip(f), 
the LHCb collaboration has recently reported Ar s = 
(0.116±0.019)ps _1 [4], which represents the current most 
precise measurement of this observable. 

As we pointed out in Ref. [5], the sizable Ar s compli- 
cates the extraction of the branching ratios of _B s -meson 
decays, leading to systematic biases as large as 0(10%) 
that depend on the dynamics of the decay at hand. 

In the case of the B® — > /i + fj," channel, the comparison 
of the experimentally measured branching ratio with the 
theoretical prediction (JlJ is also affected by this effect, 
which has so far been neglected in the literature. 



It can be included through a measurement of the effec- 
tive B° — > ^ + fJ.~ lifetime. As B® —> ^ + ^ is a rare decay, 
it turns out that this observable offers another sensitive 
probe for NP that is theoretically clean and complemen- 
tary to the branching ratio. 

II. THE GENERAL B s -+ /x+/x~ AMPLITUDES 

The general low-energy effective Hamiltonian for the 
B® — > decay can be written as 

n cS = -^y f ;y tt a[c 10 o 10 + c s o s + c P o P 

V 27T 

+ C' w O[ + C' s O' s + C' P O' P }. (2) 

Here Gp is Fermi's constant, the V qq i are elements of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix, a is 
the QED fine structure constant, the Cj, C[ are Wilson 
coefficients encoding the short-distance physics, while the 

Oio = (s lf ,P L b)(£^ l5 e) 
O s = m b (sP R b){U) (3) 
P = m b (sP R b)(£ l5 e) 

are four-fermion operators with Pl,r = (1 T 7s)/2 and 
rrib is the 6-quark mass. The 0\ are obtained from the O.j 
by making the replacements Pr, o Pr. Only operators 
resulting in non-vanishing contributions to B® — > /x + /x _ 
are included in ([2]). In particular the matrix elements of 
operators involving the l^l vector current vanish. 

This notation is similar to Ref. [6], where a model- 
independent analysis of NP effects in b — > s transitions 
was performed. In the SM, as assumed in ([!]), only C\q 
is non- vanishing and given by the real coefficient Cf^ 1 . 
An outstanding feature of B° s — > \i + [i~ is the sensitivity 
to (pseudo-) scalar lepton densities, as described by the 
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0(p)s and 0'^ s operators. Their Wilson coefficients are 
still largely unconstrained and leave ample space for NP. 

The hadronic sector of the leptonic B® — > decay 
can be expressed in terms of a single, non-perturbative 
parameter, the P s -meson decay constant Jb s OP- 

For the discussion of the observables in Section |III[ 
we go to the rest frame of the decaying meson and 
distinguish between the ^lMl anc ^ MrMr helicity config- 
urations, which we denote as Ma A*a w ith A = L, R. In 
this notation, /^Ml ano - MrMr are related to each other 
through a CP transformation: 

|(m+Ml)cp) = {CV)\nUl) = e #Cp(w) |M^R>, (4) 
where e l< ^ cp ( wi ) is convention-dependent. We then obtain 

A(B° S -> = ( f r x ^\H cS \B° s ) = -^V*V tb a 

x /B.Msm.Ce^Wf 1 -^)/ 2 [ Vx P + S] , (5) 
where Mb 3 is the B s mass, i% = +1 and ?7r = —1, and 



P = 



S 



c 



10 



m b 



2m u \ m b + m s 



M% 2m„ 



m b 



Cp — C'p 
r^sM 



(Cs-C s 



m b + mj V Cf n M 



(6) 



(7) 



The P = |P|e i¥,p and 5 = \S\e iips carry, in general, non- 
trivial CP-violating phases tpp and tpg. However, in the 
SM, we simply have P = 1 and S = (see also Ref. [5]). 
The </>cp(At/i) factor in ([5]) originates from using the op- 
erator relation (CV)^(CV) = 1 and Q in the leptonic 
parts of the four-fermion operators. 



+yxP + 5 

-T] X P* + S* 



(11) 



Here we have taken into account that the Bg-B® mixing 
phase (j> s = 2arg(V t * V t b) is cancelled by the CKM factors 
in ^ and ([9]), and that the convention-dependent phase 
^Cp[B s ) is cancelled through ([9]), whereas 0cp(mm) Slm ~ 
ply cancels in the amplitude ratio. We notice the relation 



£l& = M = 1- 



(12) 



The observables £\ contain all the information for cal- 
culating the time-dependent rate asymmetries [7]: 

T(B° s (t) -> - T(B° s (t) -> 

r(B°(t) -> Ma >a ) + r (5°(t) -> m^Ma) 

= C A cos(AM a t) + S x sin(AM a t) 
cosh(y s t/ r Bs ) + _4^ r smh(y 5 t/rB s ) ' 

Here AM S is the mass difference of the heavy and light 
B s mass eigenstates, and 



Us 



T Ba AT s /2 = 0.088 ±0.014, 



(14) 



where Tp 3 is the B s mean lifetime; the numerical value 
corresponds to the results of Ref. jj. CP asymmetries 
of this kind were considered for B s ,d — > i + £~ decays (ne- 
glecting Ar s ) in various NP scenarios in Refs. [BTHOj. 
The observables entering (fl3| are given as follows: 



C x EE 



16 1 2 

16 1 2 



2\PS\cos{(pp~^>s) 
\P\ 2 + \S\ 2 



(15) 



III. THE B B -)• fi + fi- OBSERVABLES 



_ 2Im£ A |P| 2 sin2<^-|S| 2 sin2^ s 

*A = -, , 19 — mi, , | |9 (+OJ 



1 + 16 



p 



\s\ 2 



For the observables discussed below we need the 

A(B a s ^ f ,+ f i-) = ^-^\nl ff \B s ) (8) 

amplitude. Inserting again {CVy{CV) — 1 into the ma- 
trix elements of the four-fermion operators and using 
both Q and {CV)\B° S ) = e ^ Cp ( B =) |P°), we obtain 

G 



A{B° S ^ n+ixl) = - 



x e 



i[?icp(-B 3 )+<Acp(w)(l-»7A)/2] [_ 



t/aP* + S* 



(9) 



which should be compared with ([5]). We observe that 
\MB°s ^.rMl.r)I = 1^(3? Mr.lM r . l )I- (10) 



Following the formalism to describe B®-B® mixing dis- 
cussed in Ref. [7j, we consider the observable 



6 = -e 



A(B° S 



Ma Ma. 



|X _ 2 Re £a _ |P| 2 cos2y> P - |5| 2 cos2 V? g 

|P| 2 + |S| 2 



-^Ar — 



1 + 161 



(17) 



It should be emphasized that due to |12j) Sep = S\ and 
-^Ar = -^Ar do not depend on the helicity A of the muons 
and are theoretically clean observables. 

Since it is difficult to measure the muon helicity, we 
consider the rates 

r(s°(i)^+/0= £ r(B°(t) -> /i+ M A) ! (is) 

A=L,R 

and obtain then the CP-violating rate asymmetry 

r(g°(t) -» /xVJ - r(B° a (t) -» mV-J 
r(B°(t) ^ /i+M-) + r(po(t) ^ 

S C p sin(AM s i) 



cosh(y ;j t/ t Bs ) + ^. A r smh(y s t/ t Bs ) ' 



(19) 
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where the C\ terms ( 15 1 cancel because of the f]\ factor. allows the extraction of 



It would be most interesting to measure ( 19 1 since 



a non-zero value immediately signaled CP-violating NP 
phases. Unfortunately, this is challenging in view of the 
tiny branching ratio and as tagging, distinguishing be- 
tween initially present B® and B® mesons, and time in- 



formation are required. An expression analogous to ( 19 ) 
holds also for B^ — > decays. 
In practice, the branching ratio 

1 f°° 

BR (B. -> M + M") exp ^ 2 / M B >® dt 

is the first measurement, where the "untagged" rate 
(T(B s (t) -+ /)) = T(B° s (t) -> /) + T(B° s (t) -+ f) 



-t/TB a 



[cosh(y s t/r Bs ) + Aat smh(y s t/T Bs )] (21) 



is introduced [5j |TT]. The branching ratio (20 1 is ex- 
tracted ignoring tagging and time information. As shown 
in Ref. (-1 . due to the sizable width difference, the ex- 
perimental value (20) is related to the theoretical value 



(calculated in the literature, see, e.g., Refs.p][6]) through 
1 



BR(B S -> h + li 
where 



1 + -^Ar y s 



BR(B S 



BR(B S 
BR{B S - 



± MM 
M + M~ )sm 



= \P[ 



M M )cxpj 
(22) 

(23) 



The y s terms in ( |22[ ) were so far not taken into account 
in the comparison between theory and experiment. 
^Ar depends sensitively on NP and is hence essentially 



unknown. Using (14) and varying Aat € [— 1,+1] gives 
ABR(B S -> M + /OL = ±VsBR(B s -> ^ + M~)cx P , (24) 



which has to be added to the experimental error of ( 20 ) . 

In the SM, we have A™ — +1 and rescale (nj) corre- 
spondingly by a factor of 1/(1 — y s ), which results in 

BR(B S -> m + M~)smL = (3.5 ± 0.2) x 10~ 9 , (25) 



where we have used ( 14 ). This is the SM reference for the 



comparison with the experimental branching ratio ( 20 ) 



IV. THE EFFECTIVE B s -»• fj,+ fj,' LIFETIME 



A at y s 



(1 - l£)V7 



(1 



V 2 s)tb s 



2tb s - (1 



y 2 s ) T u 



yielding 



BR{B s ^n+^-) 
BR{B S -> fi+fj,- 



(1 



/exp 



(27) 



(28) 



We emphasize that it is crucial to the above equations 
that A at in ( 17 1 indeed does not depend on the helicities 
of the muons, i.e. A at = A Ar . 

Effective lifetimes are experimentally accessible 
through the decay time distributions of the same sam- 
ples of untagged events used for the branching fraction 
measurements, as illustrated by recent measurements of 
the B° -> J/ipfo and B° -> K + K~ lifetimes [H] by 
the CDF and LHCb collaborations: both attained a 7% 
precision with approximately 500 events, while an even 
larger sample of B® —> events can be collected 

by the LHC experiments, assuming the Standard Model 
value of the P° — > /J + /i~ branching fraction. Although 
a precise estimate is beyond the scope of this article, we 
believe that the data samples that will be collected in the 
planned high-luminosity upgrades of the CMS and LHCb 
experiments [T3j can lead to a precision of 5% or better. 



CONSTRAINTS ON NEW PHYSICS 



In order to explore constraints on NP, we introduce 



BR(B S /x + M~)exp 
~ BR(B S -> (J,+n-)su 
1 + y s cos 2ipp 



1 + AArUs 



(IPI 2 + ISI 2 ) 



i-y. 2 



l-y s cos 2ip s 
i - 1/2 



\S\ 2 , (29) 



where we have used (17) and (22). Using |T]) and the 
upper bound [2] yield R < 1.4, neglecting the theoretical 
uncertainty from R). In the case of y s = 0, R fixes a 
circle in the |S plane. For non-zero y s values, R 
gives ellipses dependent on the phases (fp t s- As these 
phases are in general unknown, a value of R results in 
a circular band. We obtain the upper bounds \P\, \S\ < 
yj (1 + y s )R- As R does not allow us to separate the S 
and P contributions, there may still be a large amount 
of NP present, even if the measured branching ratio is 
close to the SM value. 

The measurement of T^+u- and the resulting observ- 
able Aaf allows us to resolve this situation, as 



With more data available, the decay time information 
can be included in the analysis. As we pointed out in 
Ref. [5], the effective lifetime 

fZ°t(T(B a (t)^ dt 
V>- = r°°/wn 77\ . ..4-..-^ ^ ( 26 ) 



f~{T(B s (t)^ dt 



\S\ 



|P| 



/ COS 2(fp — A AT 

cos 2(ps + Aaf 



(30) 



fixes a straight line through the origin in the |P|-|S| 
plane. In Fig. [TJ we show the current R constraints in 
the IPj-jS 1 ! plane, and illustrate also those corresponding 
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0.4 0.6 0.8 
\P\ 



FIG. 1: Current constraints in the | jF-'I — | jS'] plane and illus- 
tration of those following from a future measurement of the 
effective B s — > lifetime yielding the Aaf observable. 
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FIG. 2: Illustration of allowed regions in the R-A&t plane 
for scenarios with scalar or non-scalar NP contributions. 



to (30 1. In Fig. [2j we illustrate the situation in the ob- 
servable space of the R-Aat plane. It will be interesting 
to complement these model-independent considerations 
with a scan of popular specific NP models. 

Let us finally note that the formalism discussed above 
can also straightforwardly be applied to B s u) — > t + t~ 
decays where the polarizations of the r leptons can be 



inferred from their decay products |10j . This would allow 
an analysis of (13 1, where non- vanishing C\ observables 



would unambiguously signal the presence of the scalar S 
term. Unfortunately, these measurements are currently 
out of reach from the experimental point of view. 

VI. CONCLUSIONS 

The recently established width difference Ar,, implies 
that the theoretical B® — » branching ratio in Q 

has to be rescaled by 1/(1 — y s ) for the comparison with 
the experimental branching ratio, giving the SM refer- 
ence value of (3.5 ± 0.2) x 10" 9 . The possibility of NP in 
the decay introduces an additional relative uncertainty 
of ±9% originating from Aat € [—1, +1]. 

The effective B s — > /i + fi~ lifetime t„+„- offers a new 
observable. On the one hand, it allows us to take into 
account the B s width difference in the comparison be- 
tween theory and experiments. On the other hand, it 
also provides a new, theoretically clean probe of NP. In 
particular, t m + p - may reveal large NP effects, especially 
those related to (pseudo-)scalar £ + £~ densities of four- 
fermion operators originating from the physics beyond 
the SM, even in the case that the B® —> branching 
ratio is close to the SM prediction. 

The determination of t„+u- appears feasible with the 
large data samples that will be collected in the high- 
luminosity running of the LHC with upgraded experi- 
ments and should be further investigated, as this mea- 
surement would open a new era for the exploration of 



B, 



at the LHC, which may eventually allow 



the resolution of NP contributions to one of the rarest 
weak decay processes that Nature has to offer. 
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